PLANE-SYMMETRIC SPACETIMES 
WITH POSITIVE COSMOLOGICAL CONSTANT. 
THE CASE OF STIFF FLUIDS 



PHILIPPE G. LEFLOCH 1 AND SOPHONIE B. TCHAPNDA 2 

Abstract. We consider plane-symmetric spacctimes satisfying Einstein's field 
J_j 1 equations with positive cosmological constant, when the matter is a fluid whose 

pressure is equal to its mass-energy density (i.e. a so-called stiff fluid). We 
study the initial-value problem for the associated Einstein equations and es- 
tablish a global existence result. The late-time asymptotics of solutions is also 
rigorously derived, and we conclude that the spacetime approaches the de Sit- 
ter spacetimc while the matter disperses asymptotically. A technical difficulty 
dealt with here lies in the fact that solutions may contain vacuum states as 
well as velocities approaching the speed of light, both possibilities leading to 
singular behavior in the evolution equations. 

u 

1. Introduction 

The study of global properties of cosmological spacetimes is a fundamental prob- 
lem in mathematical relativity, as it provides a first step toward understanding 
fundamental issues such as the structure of singularities and the cosmic censorship 
conjecture. Such a study can be reduced to investigating the global existence and 
asymptotic behavior of solutions to the Einstein equations, possibly coupled to the 
equations of motion for a specific matter model. In the present paper, we treat 
the class of perfect fluids whose pressure p and mass-energy density fj, > coin- 
cide. This is a limiting case (7 = 2) with the class of pressure laws p — (7 — 1) fj,, 
in which 7 G [1,2] is referred as the adiabatic exponent of the fluid. Our main 
result concerns the initial- value problem for the associated Einstein equations: we 
establish a global existence result and rigorously determine the late-time asymp- 
totic behavior of solutions. This allows us to conclude that the spacetime is future 
geodesically complete and approaches the de Sitter spacetime whereas the matter 
asymptotically disperses. 

Observe first that singularities generically arise in initially smooth solutions to 
the fluid equations, that is, shock waves in the general case 7 £ (1,2] and shell- 
crossing singularities in the case 7 = 1. This is true even when gravitational effects 
are taken into account [4]. If the solution is to be continued beyond shock waves, 
it is necessary to lower the regularity of initial data and search for weak solutions, 
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as investigated by LeFloch and co-authors (cf. the review [T] and the references 
therein). 

On the other hand, existence of smooth solutions even in a long-time evolu- 
tion can sometime be established in physically interesting situations. This is es- 
pecially true when a cosmological constant is included, as we do in the present 
paper. Global-in-time solutions and the existence of future geodesically complete 
spacetimes can be established under a smallness condition on the initial data, as 
recognized by Tchapnda |10j for 7 = 1 and under the assumption of plane symme- 
try and, later, without symmetry and for 7 € (1, 4/3), by Rodnianski and Speck [5 
and Speck [7]. 

As far as the limiting case 7 = 2 is concerned, plane symmetric spacetimes 
have been investigated by Tabensky and Taub (8| and LeFloch and Stewart [2]. In 
particular, [8] relies on two different coordinate systems in their analysis, a comoving 
coordinate system in which the fluid is at rest, and a characteristic coordinate 
system. On the other hand, the work [UH] introduced the notion of weakly regular 
solutions to the Einstein equations. 

In the present paper, we rely on areal coordinates, a coordinate system in which 
the time is defined to be the area-radius function determined by surfaces of sym- 
metry. In these geometry-based coordinates, we prove a global-in-time existence 
theorem (in the future direction) for plane-symmetric solutions to the Einstein-stiff 
fluid equations with cosmological constant. Importantly, we also derive the leading 
asymptotic behavior of solutions and conclude with the future geodesic complete- 
ness of the constructed spacetime. 

Our analysis relies on a change of fluid variables that allows us to write the fluid 
equations in a way analogous to the case of a massless scalar field, and then to take 
advantage of techniques for semi-linear hyperbolic equations. (A similar structure 
was observed in [11].) A specific technical difficulty overcome in this work originates 
in the fact that solutions may naturally contain vacuum states as well as velocities 
approaching the speed of light, both possibilities leading to singular behavior in the 
evolution equations. 

Note finally that our results extend to compressible fluids the conclusions ob- 
tained by Tchapnda and Rendall [9] for the Vlasov equation of (collision-less) kinetic 
dynamics. 

The outline of the paper is as follows. Section 2 is concerned with the derivation 
of the field equations for stiff fluids under plane-symmetry. Next, in Section 3 we 
develop the local existence and uniqueness theory and then, in Section 4, determine 
the global geometry and asymptotic behavior of the spacetimes under consideration. 



2. ElNSTEIN-STIFF FLUID EQUATIONS 

Gravitational field equations. We consider spacetimes (M, g) such that the 
manifold has the topology M = I x T 3 , where / is a real interval and T 3 = 
S 1 x S 1 x S 1 is the three-torus. The metric g and the matter fields are required to 
be invariant under the action of the Euclidean group E2 on the universal cover. It 
is also required that the spacetime has an ^-invariant Cauchy surface of constant 
areal time. In such conditions the metric can be expressed in the form 



(2.1) 



ds 2 = _ e 2v(t, X ) dt 2 + e 2X(t, X ) dx 2 + t 2 {dy 2 + ^ 
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where the time variable describes t > and the spatial variable the interval x £ 
[0,1], while the variables y and z range in [0,2-7r]; the metric coefficients rj and A 
are periodic in x with period 1. The Einstein equations read 

(2.2) G^+Ag" 13 =8nT al3 , 

where G Q/3 is the Einstein tensor, T Q/3 the energy-momentum tensor and A is the 
cosmological constant which we assume to be positive. We also introduce the 
notation 

p = e 2 "T 00 , j = e x+r >T 01 , S = e 2X T 1 \ p = t 2 T 22 

which defines the fluid variables of interest. 

After a tedious computation in the above coordinates, (|2.2p take the form of 
the following evolution and constraint equations (where the subscripts i, x denote 
partial differentiation) : 

(2.3) e' 2n {2t\ t + l)-At 2 =87:t 2 p, 

(2.4) e- 2ri (2tf] t - 1) + At 2 = 8vrf 2 S, 

(2.5) Vx = -47rfe A +"j, 

(2.6) e" 2A (r] xx + 7] x (r] x - X x )) - e" 2 " \X tt + (X t - m )(X t + i)J + A = 8np. 

Stiff fluid equations. The so-called stiff fluid under consideration is an isentropic 
perfect fluid with energy density fi > equal to its pressure, that is, p = p. The 4- 
velocity vector U a of the fluid is normalized to be of unit length: U a U a = — 1. The 
plane symmetry allows us to set U a := (,(e~ 71 , e~ x u, 0, 0), where £ — (1 — u 2 ) -1 / 2 
is the relativistic factor and u is the scalar velocity satisfying |tt| < 1. The energy 
momentum tensor for the stiff fluid is 



T ap = p(2U a U p +g aP ), 



that is 



(2.7) 



\ — u l \ — u z 

T 11 = e - 2X ^-,p =: e- 2A 5, T 22 = T 33 - i" 2 ^, 

1 — w 



while, due to the above assumptions, all the other components vanish identically. 
The stiff fluid equations read 

(2.8) V a T afi = 0. 

We can assume that the components T a2 and T" 3 vanish identically, while by 
computing the remaining two components we arrive at the two evolution equations 

Pt + e v - x Jx = -2X t p - 2 Vx e r '- x j -~(fi + fi), 
(2-9) I 
Jt + e r '~ V = -2X t j - 2 Vx e r '- x p - -j. 
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The equations may be put into a simpler form, as follows. Observe that the 
first-order principal part of (|2.9|) is a strictly hyperbolic system of two equations 
associated with the two distinct speeds ±e v ~ x . Introducing the Riemann invariants 

1 + u 1 — U 

2.10) r:=- fi = p + j, s:=— — fJ , = p-j, 

1 — u 1 + u 

and the directional derivatives 

D+ := d t + e^ x d x , D- := d t - e^ X d x , 
and then combining the equations in (|2.9I) together, we obtain 

D+r = -2 [Xt + n x e^ x + f) r - \ ^JFs, 

(2.H) K y 2 

D~s = -2 f At - n x e^- x + - J s - ^V^- 

Finally, the expressions for Xt and r\ x taken from (|2.3p and (12.51) can be plugged in 
(|2.11l) . and by setting X = e n y/r and Y — e v y/s we arrive at 

D+X = -Ate 2ri X - -Y, 
(2.12) i 

£>-y = — X - Ate 2r >Y, 
t 

which we will refer to as the stiff fluid equations for the unknowns r and s. 

Basic properties. It is easily checked that (|2.6p is a consequence of the equations 
(|2.3p - (|2.5[) . (|2.9j) . One can also check that (|2.5I) is a constraint equation, that 
is, it is automatically satisfied for all times once it is satisfied on an initial Cauchy 
hypersurface. Therefore, we will work with (|2.3[) . (12.4[) and (|2.12[) for the unknowns 
rj, A, r and s. Observe that by definition r and s must be non-negative. From the 
definition we see that S = p = (r + s)/2. 

We will solve the initial-value problem with data prescribed on the hypersurface 
t = 1. Observe that once the fluid variables have been determined, the metric 
coefficient 77 is obtained by integrating (|2 .4[) in the time direction, i.e. 



(2.13) e — 2 n(t,x) _ j + _ / r 2(A-47r(r + s)(r,a;))dr 



with rj := 77(1 , •). Next, 77 being known, the following equation (obtained from (|2.3p 
and ED), 



(2.14) A t (4,1) =77 t (^)+Ate 2 "- i, 
is integrated in time to yield the second metric coefficient 

(2.15) X(t,x) = X(l,x)+ J X t (T,x)dT 

with A = A(l, •). Therefore, it will be enough to concentrate on the stiff fluid equa- 
tions (|2.12p together with the metric equation (|2.13p . that determine an evolution 
system for the unknows 77, r, s. 

Observe that there exists some T* > 1 such that the right hand side term in 
(|2.13p is positive on [1, T*) x [0, 1]. Estimates for r and s can easily be derived as 
follows. The expressions for At and 77^ taken from (|2.3p and (|2.5p can be plugged 
in (f2~TTj) to yield 
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D+r = -(87rfe 2, 's + Ate 2 '' + -) r - ~Vrs, 

(2.16) * * 
D~s = -(87rte 2 "r + Ate 2ri + -) s - -^/Ts. 

Using the fact that r and s are positive, this implies 

D + r < -t^r, D~s < -t~ X s, 

and integrating this along the characteristic curves associated with the operators 
D ± implies that 

(2.17) r<r(l,-)i - \ s< s(l,.)t -1 . 

As a consequence, if 77 is bounded then so are X and Y. 

A straightforward computation leads to the following result. 

Lemma 2.1. Set 

h = (A - v) x ^- x - Ate 2 ", b 2 = -2Atr/ x e 2 ", 

ba^-A^e^-Ate 2 *, 6 = -l. 
// A and y sofoe i2.1ty) then X x and Y x satisfy 

, x D+X x = hX x + bY x + b 2 X, 

2.18) 

D-Y x = bX x + b 3 Y x + b 2 Y. 
The following result will be used to obtain bounds on derivatives of X and Y. 

Lemma 2.2. Set 

K{t) = sup{(A + Y)(t, x)\xe [0, 1]}, 
A(t)=sup{(\X x \ + \Y x \)(t,x) | x e [0,1]}, 

v(t) = sup{|(A - fi) x \e r '- x + Ate 271 + - | x G [0, 1]}, 

/i(t) = 2Atsup{|77 K |e 2? ' | i G [0, 1]}. 

If(X,Y) and (X X ,Y X ) solve WWl and 1%J8\) , respectively, with X x (l) = (e^y/f) x 
and Y x (l) = (e v \^) x , then 

(2.19) A(t) < A(l) + J (v(t)A(t) + h(T)K(r)) dr. 

Proof. Equations (|2.18l) can be written in the form 

^-X*(t,7i(*)) = (b 1 X x + bY x + b 2 X)(t, ll (t)), 
jY x {tn2(t)) = (bX x + b 3 Y x + b 2 Y)(t, l2 (t)), 

where 71 and 72 are the integral curves corresponding to D + and D~ respectively. 
Integrating this over [l,i], taking the absolute value in each equation, adding the 
resulting inequalities and taking the supremum of each term yields (12 . 19[) . □ 
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3. Local existence theory 

Main statement of this section. We are interested in regular solutions, denned 
as follows. 

Definition 3.1. A regular solution to the plane- symmetric Einstein- stiff fluid equa- 
tions consists of two metric coefficients ry, A and Riemann invariants r, s given as 
continuously differ entiable functions defined on [1,T] x [0, 1] and periodic in space. 

We pose the initial- value problem by choosing some functions fj, A, f, s as periodic 
functions on [0, 1] satisfying the constraint 

(3.1) n x = -2Trte J+ v(r-s), 
and, on the initial hypersurface t = 1, we impose 

(3.2) (r 7 ,A,r, S )(l,-) = (r?,A ) r,s). 

Theorem 3.2 (Local existence and uniqueness theory in the Riemann variables). 

Given periodic, continuously differ entiable data rj,\,r,S prescribed on the initial 
hypersurface t = 1 and satisfying the constraint f ff. 1]) . there exists a future devel- 
opment which consists of continuously differentiable functions 77, A, r, s defined on 
some time interval [1,T) (with T € (l,oo])^ that are periodic in space and satisfy 
the stiff fluid equations iS.lty) , together with the evolution equations i2.3\) and ^-4ty - 

Once the Riemann invariants r and s are known, the primary fluid variables n 
and u can be determined from equations (|2.7[) and (|2.10l) : 

/' = 



By construction, the Riemann invariants are bounded, and this property is equiv- 
alent to the following restriction in the fluid variables: 

(3.3) i|%<l. 

1 =p it 

Observe that this condition allows the density to vanish, and the velocity component 
it to approach ±1, which is the normalized light-speed. The condition is equivalent 
to saying 

(3.4) 0<^<1-M 2 - 

Theorem 3.3 (Local existence and uniqueness theory in the fluid variables). Un- 
der the assumptions of Theorem \3.'A the problem with initial data satisfying the 
uniform bound (13.41) admits a local-in-time solution which is unique in the follow- 
ing ( generalized) sense: if [i\ , U\ and fj,2, «2 denote fluid solutions to the same initial 
value problem, then 

either fii = yU 2 > and U\ =1*2, 
or ijli = fj,2 = and U\ , u 2 are arbitrary. 

Proof of the local existence result. We rely on an iterative argument and define 
a sequence (r/ n , r n , s n ) in the following way. 

(1) For t e [1, +00) and x € [0,1], we set (rjo, ro, so)(t, x) :— (rj,r,s)(x), T — 
+00. 
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(2) If T)n-i, r n -i, s n -i are regular on [l,T n _i) x [0,1] with T n -i < oo, then 
we define T n to be supremum of all t' € (l,T„_i) such that 

e -2fj(x) 1 ft 

t 2 (A - 47r(r„_i + s n _i)(r,x)) dr > 



for all x € [0, 1] and t e [1, f'], and we then set 

p —2rj(x) i ft 

(3.5) e -a»»(t,*) = + _y r 2( A -4 7 r(r„_ 1 + Sn _r) (r,x)) dr. 

(3) We define r n and s n such that X n = e^-v/r^, and V n = e J) "y / s^" are solu- 
tions of the system 

Pn-l^n = a n-lX n -l +bY n -i, 

(3.6) 

where a n -i = — Ate 2 * 1 ™- 1 , b = — j. is the D ± -operator corresponding 
to the n-th iterate. We prescribe the same initial data (13.21) for all n. 
Observe that T n > T* for all n, so that all the iterates are well-defined and 
regular on the fixed time interval [1,T*). 

In order to prove that the sequence of iterates converges to a regular solution, we 
establish uniform bounds on the iterates as well as their time and space derivatives, 
and we prove their uniform convergence. This is done in a series of lemmas. 

In the sequel we denote by || || the sup-norm on the function space of interest, 
C denotes a constant that may change at each occurrence. 

Lemma 3.4. The sequences r\ n , X n , Y n , r n , s n and (j] n ) t are uniformly bounded 
in n, in the sup-norm by a continuous function oft, on a time interval \\,T^\. 

Proof. Set 

P n (t) := sup{e 2 ^^) \x e [0,1]}, 

K n (t) := sup{(X n + Y n )(t,x) \x€ [0,1]}. 

Using equations (|3.6p , we apply the same argument used in the proof of Lemma 
and obtain 

(3.7) K n {t) < K Q + J m n _i(r)if„_i(r) dr, 
with 

m n (t) = sup{Ate 2 ''" + -; x € [0, 1]} 
<f(l + A)(l + P„(f)), 

so that 

(3.8) K n {t) < K + (1 + A) [ r(l + P„_ 1 (r))^„_ 1 (r) dr. 



On the other hand equation (|3.5|) implies 

(3.9) ( Vn ) t = i - |te 2 "» + 27rte 2 ^- 2 ^(X 2 _ 1 + Y 2 ^), 
and since e" 2 ''— 1 < e ^+ At3 < C(l + A)t 2 , it follows that 

(3.10) P n {t) < ||e 2 l + C(l + A) f r 3 (l + X„_ 1 (r)) 2 (l + P„(r)) 2 dr. 



s 
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Now defining Q n {t) :— sup{K m (t) + P m (t); m < n} and adding (|3.7[) and (|3.10p we 
arrive at 

(3.11) Q n {t) <K + \\e^\\ + C{1 + A) J r 3 (l + Q„(r)) 4 dr. 

Let [1,TW) (with T (1) G (1,T*]) be the maximal interval of existence for the 
solution z\ of the integral equation 

z x {t) =K + \\e^\\ + C(l + A)^T 3 (1 + ^i(r)) 4 dr, Zi(l) = K + \\e^\\. 

Then Q n (t) < for all n G N and t G (1,T( X )). The same is true for K n and 

P n . It follows that ?7„, A"„, Y„, and then r„, s„ and (?? n )t are uniformly bounded. 
To bound (r? n )t, we use (|3.9I) . □ 

Lemma 3.5. TTie sequences (i] n ) x , (X n ) x , (Y n ) x , {X n ) t , (Y n ) t , {r n ) x , (s n ) x , (r„) t 
and (s n )t are uniformly bounded in n, the sup-norm by a continuous function oft 
on a time interval [1,T^]. 



Proof. Set 



A n [t) 
A 
B n (t) 



■■sup{\(X n ) x \ + \(Y n ) x \(t,x) \xG [0,1]}, 
= sup{(|X x | + |F x |)(i) \xe [0,1]}, 
: S up{\(e- 2 ^)) x \ \ x g [0,1]}. 
Then taking the spatial derivative in (|3.6|) gives the following equations: 

£>+-l(*n)* - Vn-l)^- 1 -^- 1 ^^ - 2A%„_ 1 ) x e 2r '»- 1 X tl 

£»-_ 1 (F n ) :c =(r? n _i - A^i^e"- 1 -^- 1 ^)* - 2Ai(7 ?rl _i) ;c e 2 ""- i y n 

But, using Lemma l3.4| we have 

|(A„_i -T] n -i)x(s)\ = |(A-?7)x+2A y r(77„_i) a; e 2 ' ) "- 1 dr| 
<C« 2 (1+B n _i(*)), 



so that applying the same argument as in Lemma 12.21 and using Lemma 13.41 again, 
we obtain 

(3.12) A n (t)<A Q + cj^ r 2 (l + S„_ 1 (r))(l + A„_ 1 (r)+A l (r))) dr. 
On the other hand we have 

(3.13) {e- 2 ^>*\ = ~ 2Vx t € 2 " j\ 2 (r n ^+ 3n ^) x (T,x)) dr, 
which implies 

(3.14) B n (t) < 21177^-^11 T 2 (A n -i + B n ^){r) dr. 
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We have used the fact that 

\(r n + SnU = \{e- 2Vn ) x {X 2 n + Y^) + 2e- 2 ^{X n {X n ) x + Y n {Y n ) x )\ 
< C(A n + B n )(t). 

Now defining E n {t) := sup{A m (t) + B m (t);m < n} and adding (|3~T2|) and (|3T4|) 
we arrive at 

(3.15) E n (t) <A + 2\\n x e-^\\+C^ r 2 (l + E n (r)) 2 dr. 

Let [1,T^ 2 ') (with T^ 2 ' < T^) be the maximal interval of existence for the solution 
Zi of the integral equation 

z 2 (t) = A + 2||^e- 27T || + cj r 2 (l + z 2 (t)) 2 dr, 

z 2 (l) = A + 2\\rj x e- 2 % 

Then E n (t) < z 2 (t), for all n e N and t e (1,T (2) ). The same is true for A n and 
B n . It follows that (r) n ) x , (X n ) x , (Y n ) x , (X n ) t , (Y n ) t , (r n ) x , (s n ) x and then (r n ) t , 
(s n )t and {rj n ) tx are uniformly bounded. □ 

Lemma 3.6. The sequences irj n ), (X n ), and (Y n ) converge uniformly on [1,T^] 
for allTW less than T {2 \ 

Proof. For t € [1,T^], define 

6 n (t) :=sup{\X n+ i-X n \(t,x) + \Y n+1 -Y n \(t,x);xe [0,1]}, 

a n (t) :=sup{\\(r] n+1 - r] n )(s)\\ + \\(X n+1 - X n )(s)\\ + \\(Y n+1 - Y n )(s)\\;s e [l,t]}, 

X n -=X n+ i — X n , Y n := Y n+ i — Y n . 
Combining equations (|3.6j) written for n + 1 and n gives 



(3.16) 
with 





Dn-Xn — QfiXn- 




+ i r n, 






D-Y n = bX n _x 


+ a n Y n -i 


+■ G„, 




F n = -(e 2 "" 




- ( e ""- A " 




- 1 )(^n)x, 


G n = -(e 2 "" 




f (e""- A " - 


_ e 77»-l-An 


- 1 )(* r n)«- 



Reasoning as in the proof of Lemma 12.21 we have 

(3.17) 6{t) < J (m„(T)0 n _i +sup{|F„(r,x)| + \G n (T,x)\;x € [0,1]}) dr, 
and this implies that 

(3.18) \X n \ + \Y n \ < C J a„_!(r) dr, 

we have used the mean value theorem to handle the terms e 2l) ™ — e 2r,n ~ 1 and 
e %>-A„ _ e r)„_i-A„_i^ anc j p rev i ous lemmas. 

On the other hand equation (|3 . 9[) implies 

( 7?n+1 - , M)t = - |t(e 2 ^ - e 2 "") + 27rte 2 "» + - 2 ^ ((X 2 +1 - X 2 ) + (Y 2 +1 - Y 2 )) 
+ 2nt{e 2 ^ +1 - 2 ^ - e 2r '"- 2r '- 1 )(X 2 + Y 2 ), 
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and using Lemma 13.41 and the mean value theorem it follows after integration in 
time that 

\r)n+i - Vn\ < C J (|?7 n+ i - r] n \ + \rj n - r]„_i| + \X n+1 - X n \ + \Y n+1 - F„|)(r) dr, 
so that 

(3.19) \Vn+l-Vn\ <C / ( )(r) dr. 
Combining ([3~T8]) and (|3~19| leads to 

a n (t) < C (a n + a n -i)(r) dr, 
which, by Gronwall's inequality, implies 

a n (t) < C J a„_i(r) dr, 

and by induction 

C n+1 

Oi n {t) < j— , 

n! 

and so a n — > as n — > oo. This establishes the uniform convergence of rj n , X n , and 
Y n . □ 

It follows from (13.91) that the sequence (rj n )t converges uniformly as well. In the 
following lemma, the uniform convergence of other iterates derivatives is proven. 

Lemma 3.7. The sequences (r] n ) x , (X n ) x , (Y n ) x , (X n ) t and (Y n ) t converge uni- 
formly on [1,T(% where [1,T^] C [1,T^]. 

Proof. We set 

(3 n (t) :=sup{||(r7„ +1 -r ? „) :r (s)|| + ||(X ra+1 -X„) x (s)|| + ||(y n+1 -y n ) x (s)||;s e [l,t]}. 
Taking the space derivative in equations (|3.6I) gives 

(3.20) D+(X n+1 ) x = C n , D-(Y n+1 ) x =D n 
with 

C n = (A„ - Vn) x e^- x -{X n+1 ) x - 2At(r ]n ) x e 2 ^X n - Ate 2 ^(X n ) x - j(Y n ) x 

D n = (vn - K) x e^- x -{Y n+1 ) x - 2At(n n ) x e 2 ^Y n - Ate 2 ^(Y n ) x - ^(X n ) x . 

Let 7^ and 7 2 be the integral curves corresponding to and D~ respectively, 
that start from the point (s, x) that is, for each n, 

(3.21) ( 7 i) t = e ^~ A ", ( 7 , 2 J t = - e ^~ A ", 7 i( S ) = 7 2 ( S ) = 

Integrating the first equation in (I3.20j) along 7 ^, the second one along j 2 yields 
after subtraction 

(X n+1 - X n )(s) = [ (C„(t, 7 ^(t)) - C„-i(t,7£_i(t))) dr, 

(3.22) Jl s 

(Y n+1 - Y n )(s) = / (D n (rnUr)) - A.-ifoTi-iM)) 
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But we have 

< \C n (rnUr)) - ^(t.tLiWI + \(C n - C n ^)(r)\. 
Given now any e > 0, we find, for any sufficiently large n, 
(3-24) \C n (r,^(r)) - C„(t, 7 ^i(t)| < Ce, 

we have used the uniform convergence of rj n , the uniform continuity of C n over the 
compact set [1,T^] x (^([1, T^]) U 7^_i([l, T' 4 '])), and the following inequality 
which follows from (|3.2i| 

(3.25) ^-Ti^lM^CsupjlKe^-e 2 ""- 1 )^)!! ; te[l,T^}}. 

For the second term of the right hand side in (|3.23|) we have 

C n — C n _i 

= ((A n - ij n ) x - (A„_i - r)n-i)x)e nn ~ Xn (X n+1 ) x 

+ (A„_! - TM-O^e^-^C^+j - + (e""- A " - e^- x ^){X n+1 ) x ) 

- 2At( Vn - Vn -i) x e 2 ^X n - 2M{f ]n ^ 1 ) x (e 2ri -{X n - X n _ x ) + (e 2 "» - e 2 ^){X n . 

- A<(e 2j) " - e 2 ^)(X n ) x - Ate 2 "- 1 ((X n ) x - (JC n _ 1 ) x ) - i((y„) a - (Yn-i)*), 
and 

(A„ - ry„) K = (A - 77)2; + 2A / t(-q n ) x e 2nn dr, 



so that 

|(A„ - - (A„_! - rjn-O^I < Ce + Csup{||(r;„ - r} n -i) x (t)\\ € [l.T^]}. 
Thus, for n sufficiently large, 

(3.26) ||(C'„-C , „_ 1 )(t)|| ^Ce + C^ + ^OM. 

It then follows from ((3~22l) - ((3?24l) and ([3T26)) that for n sufficiently large, 



|(X n+1 - X n ) x \{s) <Ce + C (/3 n + /3„_i)(r) dr, 
(3-27) "J, 1 

|(Y n+1 - Y n ) x \(s) <Ce + C (/?„ + /3„_i)(r) dr. 



On the other hand, taking the spatial derivative in (|3.9[) . subtracting the resulting 
equations written for n + 1 and n gives 

(?7n+l - Vrdtx 

= - M(r, n+1 - Vn ) x e 2 ^ - Ai^^e 2 "^ 1 - e 2 "») 

+ 4irt( Vn+1 - Ty^e 2 ^ 1 ^^ 2 + Y 2 ) - 4 7 rt(r 7n - Ty^^e 2 ^-""- 1 ^ 2 ^ + F t 2 _ 
+ 47rte 2 (""+ 1 -' ? ")((X„ - X n ^) x X n + (Y n - Y n ^) x Y n ) 
+ 47rte 2 ^-^((X n ^) x X n + (Y n ^) x Y n ) 

- 47rf e 2 (""-""- 1 > ((Jf n _ 1 ) s X n _i + (y„_i) a y„-i), 
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from this and the previous lemmas, it follows that, for n sufficiently large, 
(3.28) \( Vn+1 - VnU(s) <Cs + C [ {fi n + p n -i){r) dr. 



Combining p.27|) and (|3.28|1 . and taking the supremum over s <E [l,t] yields, for n 
sufficiently large, 

Pn(t) <Ce + C ^ + Pn-l)(T) dr, 

and by Gronwall's lemma it follows that, for n sufficiently large and t € [1,T^], 

5 n (t) < Ce, 

where 5 n (t) := sup{/3 m ,m < n}. The uniform convergence of (r] n ) x , (X n ) x , (Y n ) x , 
(X n ) t and (Y n ) t follows. □ 

Lemmas 13.61 and 13.71 allow us to pass to the limit in (|3.5p and (|3.6[) and obtain 
a regular solution (rj, X, Y) to our system on a time interval \1,T). It is easily 
checked that this solution is unique. Namely, let (rji, Xi, Yi), i = 1, 2, be two regular 
solutions of the Cauchy problem for the same initial data (rj, X,Y) at t = 1. Using 
the same argument as in the proof of iterates convergence leads to 

a{t) <C J a(r) dr, 

where a(t) = sup{\\(m - r) 2 )(s)\\ + ||(Xi - X 2 )(s)\\ + ^ -Y 2 )(s)\\ ;s e [l,t]}. 
It follows that a(t) = 0, for t € [1,T) i.e. the solution is unique. 

We have thus established the existence of a unique, local-in-time regular solu- 
tion {rj, A, r, s) to the Cauchy problem for the plane symmetric Einstein-stiff fluid 
equations written in areal coordinates. 

4. Global existence theory and asymptotics 

Global existence. We are now in a position to establish the following main result, 
which takes advantage of our assumption A > 0. 

Theorem 4.1 (Global existence theory and asymptotics). Under the assumptions 
in Theorem \3. 6 A the solution constructed therein is defined up to T = +oo, the 
spacetime is future geodesically complete, and the following asymptotic properties 
hold at late times: 

f] = -\nt{l + 0((\nt)- 1 )), X = \nt{l + 0({\nt)- 1 )), 

r^oit- 1 ), s = 0(r 1 ), 

(4-1) 1 1 

Vt = — t (l + 0(t- 1 )), A^-fl + Or 1 )), 

Vx = 0(l). 

Consequently, the generalized Kasner exponents associated with this spacetime (cf. (|4.11[) , 
below) tend to 1/3: 

Hm = i im = i im = i, 

t->oo k(£,x) t-s-oo K(t,X) t->oo K{t,X) 3 

where n — k\ denotes the trace of the second fundamental form k\ . 
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In particular, this shows that the spacetime approaches the de Sitter spacetime 
asymptotically. To establish this global result, we begin with a continuation crite- 
rion, based on the same notation as in the previous section. 

Lemma 4.2. Let [1,T) be the maximal interval of existence of solutions to the 
system under consideration. If sup{\r](t 1 x)\ \x E [0,1], t € [1,T)} < +00 then 
T = +00. 

Proof. It suffices to prove that under the assumption that 77 is bounded on [1,T), 
the same is true for r) x , rjt, X, Y, X x , Y x , X t , and Y t . First of all, by definition 
we have X = e v ^/r and Y = e v y/s and it follows from the decay inequalities 
(|2.17[) that X and Y are bounded. Next, recalling that r/ x = — 27rfe' ,+A (r — s) and 
Vt = + 2-Kt(X 2 + Y 2 ) — ^e 2?? , we find that r\ x and 77* are bounded as well. Here, 
we have used the fact that 

f* 



X(t,x) = (A- rj){x) +T}(t,x) - hit + A J Te An (T,x) dr. 

Taking the spatial derivative in this equation implies 

(\-r)) x (t,x) = (X-7})(x)+2aJ T( Vx e 2 ^)(T,x) dr, 

so that v(t), defined in Lemma |2"T21 is bounded. Rewriting (|2.19l) 

A(t) < A(l) + J (v(t)A(t) + H{t)K(t)) dr, 

and using the fact that h and K are bounded, Gronwall's lemma allows us to 
conclude that A, and then X x and Y x are bounded. Bounds on X t and Y t then 
follow from ([2~T2]) . □ 

We now prove that rj is bounded in order to conclude that T = +00. 

Lemma 4.3. The function rj satisfies 

sup {Mi, x) I /x€ [0,1], te [1,T)} <+oo. 

Proof. We can deduce from (l2~13l) that e - 2r >^ < ±^-l+ML^ i. e . 

( ' e -c + kt^ 

which provides a (negative, say) lower bound on rj. Now, let us prove that 



(4.3) 



f (e^ +x p){t,x) dx<Cr\ t€[l,T), ie[0,l], 
Jo 



which will eventually lead us to an upper bound for 77. 

Using the equations (|2.3[) . (|2.4p . and (|2.9[) . after some computations we find 

±^j\e*>+ x p)(t,x) dx^j = j\v + *p(-~-Me 2r >) dx - £ e 2r > (j x + 2 Vx j) dx 

1 2 

-e n+A /j dx. 
t ^ 

Since p > and 

1 e 2 "^ + 2 Vx j) dx = [ (e 2r >j) x dx = 0, 
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it follows that 



(4.4) j t ^ j\e^ +x p)(t, x) dx^j <\ f Q eV+X p( ~ 1 - AiV 7 ) dx. 



Thanks to (|4.2p . we have 



, 2 2t , -At 3 n 9C , 

-AtV < t < -3 + -^i" 3 , 

-c + 4* 3 ~ A 



so that (|4.4|) implies 



d 



< 4 / (e" +A p)(*,aO dx ] =4i 3 f\e ri+x p)(t,x) dx + t 4 d 



dt 



{e^ +x p){t,x) dx 



9C 



< 4i 3 y (e" +A p)(i, s) dx + t 3 J e ,1+x p( - 4 + ^ ) dx 



from which we deduce (|4.3I) by integration. 

We are now in a position to make use of the integral estimate (14.31) . Recalling 
that r] x — —4irte v+x j and < j < p, we control the spatial oscillation of r\ at each 
time, as follows: 

/>1 />1 pX />1 />1 

rj(t,x) - / rj(t,T) dr = r) x (t,a) da dr < \r) x (t,a)\ da dr 

Jo Jo Jt Jo Jo 

< Airt f (e rl+x j)(t,a) da < 4irt f (e ,1+x p)(t,a) da, 
Jo Jo 

and thanks to (|4.3[) . this implies that 



Tj(t,x)- v(t,T)dr <ct~ 3 , te [1,T), xE [0,1]. 



(4.5) 

We will have the desired upper bound on rj, provided we can control its integral. 
Recalling that r\ t — Xt = 7 — Ate 2r > and using (|4.2[) gives 



0* 



1 



At 2 



and, after integration, 



(4.6) 



2 "- A < C ■ 



< Ct~ 2 , t e [1,T), a; <= [0,1] 



Next, using ([13]), (|4T5]| . (|Q]| , and (|4T5]) . we have 



r)(t,x) dx — j rj(x) dx + J J r) t (s,x) dxds 
t ! ,1 



< C 



1 2s 7 



(l + e 2, '(87rs 2 p- As 2 )) deeds 



<C r +-lnt + 47r 



jn-\^+\ p dxds 



1 Jo 



1 JO 



— se 2ri dxds, 
2 
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thus 



1 f* 1 f* As 2 

rj(t, x) dx < C + - hit + C / s~ 5 da- - I - — r— ds 

2 Ji Mi C + ^s-s 

1 ( At \ 
< C + - In 



2 yc+±t 3 )' 

It then follows from (|4.5[) that 

-3x 1- ^ Ai 



r?(t,x) <C(l+i- 5 ) + -ln 



2 ^c + f^y' 

which leads to an upper bound for 77, i.e. 

e 2 V (t,x) < c , r a j te [ ljT ) ; [ 0jl j 5 

and the proof is complete. □ 

Late-time asymptotics. We determine now the explicit leading asymptotic be- 
havior of r, s, 77, A, At, 774 and rj Xl and then check that each of the generalized 
Kasner exponents tends to 1/3. We have proven that (see equation (12.171) ) 

(4.7) r = o{r 1 ), s = o{r 1 ). 

and the equation (|2.13l) implies 

(te- 2r >) t = At 2 ~4irt 2 (r + ,s). 
Integrating over [l,t] and using (|4.7|) . we obtain 

< Ct 2 , 

that is, e~ 2 '' = (A/3)i 2 (l + 0(i -1 )), so that 

"J 



te- 2 " - 



A t- i (i + or 1 )). 

In view of r) t = (l/2t) - (A/2)te 2r > + 27rfe 2r '(r + s), one has 

(4.8) , ft = -i(l + 0(t- 1 )) ) 

and, after integration over ?7 = — lni(l + 0((lnt) -1 )). 
Since At = T)t + Ate 2r > — (1/t), one also has 

(4.9) A t = i(l + 0(t- 1 )), 

and integrating over[l,i] gives A = lni(l + 0((ln i) -1 )) . This implies e A = O(i), 
and recalling that r\ x = — 2nte x+n (r — s) one deduces that 

(4.10) Vx = 0(1). 

Consider the generalized Kasner exponents which take the following form for the 
metric under consideration (see for instance [3]): 

n\{t,x) _ t\ t n\{t,x) _ «§(*, x) _ 1 



k(£,x) tA t + 2' K(t,x) n(t,x) tXt + 2' 
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where n(t,x) — K,\(t, x) is the trace of the second fundamental form Kij(t,x) of the 
metric. It follows from (|4.9[) that as t tends to oo, each of these quantities tends to 
1/3, uniformly in x. 

Future geodesic completeness. The late-time asymptotic expansion above al- 
lows us to establish that the spacetime is future geodesically complete, as follows. 
Let r i — y (7 q )(t) (with t = 7°(r)) be a future directed causal geodesic defined on 
an interval [1,7+) with r + maximal, and normalized so that 7°(to) = £(to) = 1 for 
some t q € [1, r + ). We are going to prove that r + = +oo. 
Since 7 is causal and future directed, we have 

g a ^l P T = -m 2 , 7t ° > 0, 

where m = if 7 is null, and m ^ if 7 is timelike. Since 4£ = 7° > 0, the geodesic 
can be parametrized by the coordinate time t. With respect to this coordinate time 
the geodesic exists on the whole interval [l,+oo) since on each bounded interval 
of t the Christoffel symbols are bounded and the right-hand sides of the geodesic 
equation (written in coordinate time) are linearly bounded in 7*, 7 2 , 7^. 
Along the geodesic we define 



w:=eH, F:=t 4 (( 7T 2 ) 2 + (7 r 2 ) 3 ). 



Using the geodesic equation it is easily checked that 

dw \ 2n-\ 1 0\2 dF 

dr dr 
The relation between coordinate time and proper time is then given by 

(4-12) % = ( 7r °- 1 - 



dt KITI ^m 2 + w 2 + F/t 2 

We will now exhibit a lower bound for dr j dt by a function with divergent integral 
on [1, +00) and, to this end, an estimate on w as a function of the coordinate time 
is needed. 

Assume that w(t) > for some t > 1. Then, as long as w(s) > 0, we have 

— = — X t w — e v ~ x r) x \J m 2 + w 2 + F/s 2 
ds 

(4.13) = 4irse 2v (j ^Jm 2 + w 2 + F/s 2 - pw) + -^w - —se 2r, w. 

Using the elementary inequality \/a + b < yfa + \fb and the equation (|4.2|) , we 
obtain 

d w 1 — As 2 e 2,? 

— < Airse 271 {\j\ - p)w H w + Airse 2v \j\^m 2 + F/s 2 . 

ds 2s 

We can drop the first two terms which are negative since \j\ < p and 

C 

1 - As 2 e 2n < ~£ s ~ 3 - 2 < 0, s sufficientlylarge, 

and we estimate the third term by Cs~ 2 (since < Cs _1 and e 2 '' < Cs~ 2 ). It 
then follows that 

(4.14) ^ < Cs- 2 . 

ds 
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Let to € [l,t) be the smallest time such that w(s) > for all s € [to,*)- Then 
integrating ()4.14|) over [io,i] gives 

w(t) < C. 

For the case w(t) < 0, it follows from (I4.13|) that, as long as w(s) < 

dw . , on/ r~5 5 TT7~o s 1 - As 2 e 2l) 

— > 47rse '( — p\/m z + + Ws z — pit) J -| tu 

as ' 2s 



> -A:irse 2ri pyjm 2 + F/s 2 + 8Trse 2ri pw 

> Cs- 2 (-l + w), 

we have used the fact that \j\ < p, 1 ~ A 2 s g e - < for large s and the elementary 
inequality y/a + b < y^a + v6- Therefore we have 

1 d(l — w) ^9 

(4.15 y — < Cs- 2 . 

1 — w as 

Let ti e [1,£) be the smallest time such that w(s) < for all s € [ti,t)- Then 
integrating (|4.15[) over [t 1: £] implies 

-«>(t) < C*. 

In either case, we arrive at 

K*)|<C, t>l. 
On the other hand equation (|4.2I) implies that 

e'' > C* -1 , f > 1, 
so we then deduce from (14.12[) that 

dr Ct- 1 

~dt - v^+WTf : 

and since the integral of the right-hand side over [1, +oo) diverges, it follows that 
r + = +oo and the proof of future geodesic completeness is completed. 
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